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lished. The three parameters, gi, 92 and gs, characterize the cross-terms in the Lagrangian, between 
the Maxwell field and terms linear in the Ricci scalar, Ricci tensor, and Riemann tensor, respec- 
tively. Static spherically symmetric equations are set up, and the three parameters are interrelated 
and chosen so that effectively the system reduces to a one parameter only, q. Specific black hole and 
other type of one-parameter solutions are studied. First, as a preparation, the Reissner-Nordstrom 
solution, with qi — q2 ~ qs ~ 0, is displayed. Then, we seek for solutions in which the electric field 
is regular everywhere as well as asymptotically Coulombian, and the metric potentials are regular 
at the center as well as asymptotically flat. In this context, the one-parameter model with qi = —q, 
q2 — 2q, q-i = —q, called the Gauss-Bonnet model, is analyzed in detail. The study is done through 
the solution of the Abel equation (the key equation), and the dynamical system associated with the 
model. There is extra focus on an exact solution of the model and its critical properties. Finally, an 
exactly integrable one-parameter model, with q\ = —q, q2 — q, q-j. = 0, is considered also in detail. 
A special sub- model, in which the Fibonacci number appears naturally, of this one-parameter model 
is shown, and the corresponding exact solution is presented. Interestingly enough, it is a soliton of 
the theory, the Fibonacci soliton, without horizons and with a mild conical singularity at the center. 
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I. INTRODUCTION 

A natural and very general extension of the Einstein-Maxwell Lagrangian yielding a general system of equations for 
a non-minimal coupling between the gravitational and electromagnetic fields, with non-linear terms, was set up and 
studied in Within this general theory, a special theory, worth of discussion, arises when one restricts the general 
Lagrangian to a Lagrangian that is Einstein-Hilbert in the gravity term, quadratic in the Maxwell tensor, and the 
coupHngs between the electromagnetism and the metric are linear in the curvature terms. The motivations for setting 
up such a theory are phenomenological, see, e.g., 0, 0] for reviews and references. This theory has three coupHng 
constants qi , (72 and 53 , which characterize the cross-terms in the Lagrangian between the Maxwell field Fij and terms 
linear in the Ricci scalar R, Ricci tensor Rik , and Riemann tensor Rikmn , respectively. The coupling constants qi , 52 
and (73 have units of area, and are a priori free parameters, which can acquire specific values in certain effective field 
theories. More specifically, the action functional of the non-minimal theory linear in the curvature is 

S^J d^x^gC, (1) 

where g is the determinant of the spacetime metric gik , and the Lagrangian of the theory is 

^ " i (f + ^ ^X^'""i^»fc^^mn) • (2) 

Here k = 2G, G being the gravitational constant and we are putting the velocity of light c equal to one, Fik = 
diAk — dkAi is the Maxwell tensor, with being the electromagnetic vector potential, and latin indexes are spacetime 
indexes, running from to 3. The tensor is the non-minimal susceptibility tensor given by 

ikmn 

qiR, 

Am ^kn An ^km \ i f Tjim ^kn jDin ^km i jDkn Am rykrn An\ , „ riikmn ( n\ 

X = ~^\9 9 ~9 9 >^~^^R 9 -R 9 +R 9 -R 9 , (3j 

where qi, 92, and (73 are the mentioned phenomenological parameters. The action and Lagrangian ([I])-® describe 
thus a three-parameter class of models, non-minimally coupled, and linear in the curvature [3, S|- Lagrangians of 
this type have been used and studied by several authors. 

The first and important example of a calculation of the three couplings was based on one-loop corrections to 
quantum electrodynamics in curved spacetime, a direct and non-phenomenological approach considered by Drummond 
and Hathrell [3]. This model is effectively one-parameter since the coupling constants are connected by the relations 
qi = —5(7, (72 — 13(7, 93 = —2(7. The positive parameter q appears naturally in the theory, and is constructed by 

using the fine structure constant a, and the Compton wavelength of the electron Ac, (7 = In these models it is 

useful to define a radius r^, an effective radius related to the non-minimal interaction, through = •\/2|g|. Thus, the 
corresponding effective radius for the non-minimal interaction in this case, is the Drummond-Hathrell radius r^^jj, 
given by t'^^jj = Ac-y/g^. In Q one also finds a quantum electrodynamics motivation for the use of generalized 
Einstein-Maxwell equations. 

Phenomenological models, i.e., models based on external considerations to obtain the coupHngs, or parameters, qi, 
g2, and (73, have also been considered. Prasanna 0,01 wanting to understand how the strong equivalence principle 
can be weakly violated in the context of a non-minimal modification of Maxwell electrodynamics, has shown that 
(7i = (72 = 0, (73 = —q, q a free parameter, is a good phenomenologically model. Another type of requirement, one with 
mathematical and physical motivations, is to impose that the differential equations forming the non- minimal Einstein- 
Maxwell system are of second order (see, e.g., jl,]^). For instance, in by imposing a Kaluza-Klein reduction to four 
dimensions from a Gauss-Bonnet model in five dimensions, thus guaranteeing second order equations for the electric 
field potential Ai, and metric gik, it was discussed a model in which (71 -I- (72 + 93 = and 2qi+q2 — 0, i.e., with qi = —q, 
(72 = 2(7 and (73 — —q. So the extra non-minimal term is a kind of Gauss-Bonnet term, and the model is called the 
Gauss-Bonnet model. Yet another type of requirement, this time purely mathematical, was suggested in [ij. The idea 
is connected with the symmetries of the non-minimal susceptibility tensor ;^'*='"" (see Eq. ([3])). For instance, one can 
recover the relations qi + q2 + qs — and 2qi + (72 = 0, used in Q, by the ansatz that the non- minimal susceptibility 
tensor Xifcmn is proportional to the double dual Riemann tensor *i?*fe„„, i.e., Xikmn = 7 *Rlkmn^ fo^" some 7 (see, [l| 
for details and motivations). Analogously, one can use the Weyl tensor Cikmn in the relation Xikmn = tuCikmn, for 
some uj, or the difference Rikmn — Cikmn instead of *Rlkmni introduce some new linear relations between (71, (72, 93, 
namely iqi -I- 92 = and 92 + 93 = 0. Yet another type of requirement is to choose the parameters so that one obtains 
exact solutions. As we will see this will lead to a model with 91 -|- 92 + 93 = and (73 = 0, i.e., 91 = —q, (72 = q, qa = 0. 
Since this model is integrable we call it the integrable model. A subcase of this has additional interest and is called 
the Fibonacci soliton. 
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Up to now we have a theory defined through Eqs. ([T])-(l3l), with each chosen set of values for the parameters qi, 
q2, and qs, giving a model. We have seen that the reduction from three-parameter models to one-parameter models, 
specified by the one parameter q and the relations between qi, (72 and 53, happens in several instances, either through 
direct calculation, as in 0, 0|, or through phenomenological and other considerations, as in [sl-fo] or and here. 
This certainly simplifies the analysis, and we will consider this one-parameter type of models, in which qi, (72, a nd 
ga, have a specified relation to the parameter q. For all these models, one can pick an effective radius = \/2|g|, as 
in the Drummond-Hathrell case, which gives the range of the non-minimal interaction between the gravitational and 
electric fields. Of course, Vq can be set to zero, in the case the world is pure Einstein-Maxwell, or otherwise can have 
a given specified value. The radius f'g^jj defined above is a candidate but in principle not the unique choice. Thus, 
possible estimations of the parameter q, and so of r^, from, for instance, astrophysical observations, are undoubtedly 
of interest (see, e.g., 

Now, after choosing a model, specified by q and by the relations between qi, 92, and (73, it is important to study 
exact solutions. Exact solutions of the equations of non-minimal electrodynamics in non-linear gravitational wave 
backgrounds were obtained in and a non- minimal Bianchi-I cosmological solution was discussed in in 

this context. Here we want to study charged black hole and other charged solutions of non-minimal models. The 
Reissner-Nordstrom solution is a standard solution in pure Einstein-Maxwell theory, with two horizons, an event and 
a Cauchy horizon, and a timelike singularity at the center (see, e.g., [lB|)- Paradigmatic charged black hole solutions 
also appear in the framework of Einstein theory minimally coupled to non-linear electromagnetic fields, as well as 
other matter fields. Such solutions were found by Bardeen and others [l3|-[2H and the main feature is that they are 
regular, without singularities inside the horizon. Within quartic gravity non-singular charged black hole solutions have 
also been found [131 ■ Since the non- minimal theory we are considering possesses new degrees of freedom, namely, the 
phenomenological parameter q and its relations to 91, (72, and (73, we beheve that these allow to introduce new aspects 
to the problem of finding black hole and other solutions of each chosen model. Three aspects can be mentioned. 
First, one wants to have a gauge in order to compare the new solutions. Thus, we study the Reissner-Nordstrom 
solution, the trivial solution in this context, where (71 = (72 = <73 = <7 with 9 = 0, in order to understand the novel 
features, such as causal and singularity structure, of the new solutions. Second, one should try to search for models 
exactly or quasi-exactly soluble. This requirement will take us, among the models cited above and the many other 
possible models, to two interesting models. They are, the Gauss-Bonnet type model where qi = ~q, (72 = and 
53 = — <7, and the integrable model where qi = —q, (72 — q, qs — 0. In both models we perform a detailed analysis. 
Third, we vary q within each of the two nontrivial models and in trying to consider non-minimal extensions of the 
Reissner-Nordstrom solution, we search for features that are similar or distinct from the two paradigmatic solutions, 
the Reissner-Nordstrom solution itself and the Bardeen solutions. In our search for charged black hole solutions in 
non-minimal models we work in Schwarzschild coordinates and impose certain requirements. The first requirement 
is connected with the electric field E{r). We demand it is a regular function on the interval < r < 00, the origin 
being also regular (i.e., the value E(0) is finite), and for large values of r the electric field is Coulombian, E{r) 
where Q is the electric charge. The second requirement is concerned with the metric functions gik- These should take 
finite values at the center {gik{0) 7^ 00). Horizons at r 7^ are not excluded, and far away the solutions should be 
asymptotically fiat. Upon these conditions and solving the non- minimal equations in certain cases we will find electric 
charged solutions with one horizon only, thus causally distinct from the Reissner-Nordstrom. However, like Reissner- 
Nordstrom, the solutions have a singularity at the center, although here the singularity is spacelike instead, as the 
Schwarzschild case, but conical, thus much milder. We have also found in one model a gravitational charged soliton, 
without horizons, where the fields are well behaved, apart from a mild conical singularity at the center. Although this 
and other solutions with horizons are almost regular at the center, we have not obtained strictly non-singular black 
hole of the type found by Bardeen and others [l3|-j23|- The difference is based on two aspects. Firstly, we consider 
spherically symmetric static solutions with goo{T)grr{r) 7^ 1, in contrast to the Schwarzschild, Reissner-Nordstrom, 
and minimal regular Bardeen solutions. Secondly, we assume that the values (?oo(0) and grr{0) are finite, but can differ 
from one. Moreover, we admit, that goo(O) can, in principle, be equal to zero. This means that the scalar curvature 
invariants for such a metric can take infinite values at the center, and the solution of the non-minimal Einstein 
equations is not regular at the center in this general sense. However, in many cases the singularity is a conical one, 
so much milder than the nasty ones of Schwarzschild and Reissner-Nordstrom. In summary, we find charged black 
hole solutions different in horizon structure from the Reissner-Nordstrom but similar to Bardeen black holes, and 
although singular, in certain cases can be considered quasi-regular conical singularities, in-between Schwarzschild and 
Reissner-Nordstrom types of singularities and the no singularities of Bardeen. 

This paper is organized as follows. In Section UIl in particular in subsection III A\ using the Lagrangian formaHsm 
of the Introduction, we establish a three-parameter non-minimal Einstein-Maxwell model. In subection III B| we set 
up static equations for studying black holes and reduce the three-parameter model to a one-parameter model. In 
Section [Ull we study specific spherically symmetric one-parameter solutions. In subsection IIII Al we define the basic 
quantitities and basic variables. In subsection IIIIBI we display the Reissner-Nordstrom solution as a preparation. 



4 



where qi = q2 = qs = q with q = 0. In subsection IIII CI we analyze in detail a one-parameter model, the Gauss-Bonnet 
model, with qi = ~q, q^ = 2q, q^ = —q (i.e., qi + q^ + qz =0 and 2qi + q^ = 0), using the known solution of the 
Abel equation, the key equation of the model, and the dynamical system associated with this model. We display the 
charged black hole solutions, and we focus on a specific exact solution and its critical properties. In subsection IIIIDI 
we consider in detail an exactly integrable one-parameter model, the integrable model, with qi = —q, q2 — q, qa — 
(i.e., 91+92+93 = and 93 — 0). We examine a special sub- model, the Fibonacci soliton, of this one-parameter 
model and present the corresponding exact solution. In subsection IIII El we present, by means of a table, the summary 
of the results of the models studied. In Section ITVl we conclude. 

II. NON-MINIMAL COUPLING, LINEAR IN THE CURVATURE, BETWEEN GRAVITY AND 
ELECTROMAGNETISM: EQUATIONS AND REDUCTION FROM THREE PARAMETERS TO ONE 
PARAMETER FOR STATIC SPHERICALLY SYMMETRIC SYSTEMS 

A. The three-parameter model: general equations 

The variation of the Lagrangian ([I])-© with respect to the metric yields (see, jl| for details), 

fTf^+'?ii;^^+rf+'?3lf^l- (4) 



Rik — -^R 9ik — K 



The energy-momentum tensor of the pure electromagnetic field, T^i!\ is 

rpia) _ 1 p prnn _ p p m /K^ 

The definitions for the other three parts of the stress-energy tensor, , T-^^^ and tI^\ are 

T.^k^ - R T^k^ - \R^kFmnF^'' - ^ff.fc V V, (J^„„F™") + ^ V, Vfc (F„,„F™") , (6) 

T^^^ = -\g^k [V„,V,(F'""i^'„) - i?,™i^""i^'„] - F'^{RaFkn + RkiF^n) ~ 

-R^''F,^Fkn-\V^i{F,nFk^)+]^'^i [V,(f^fc„F'")+Vfc(F„i^'")] , (7) 

^ik — ^dikii rmn-ris-^-r (-t'l riknls+-rk ninls)~-j^y n{i' i ^k +^k )■ \P) 

In addition, the non- minimal electrodynamics associated with the Lagrangian l[I|)-ll3]) obeys the equation 

Vfc [F''^ + x^'=™"J^™„] = , Vfei^"'= = , (9) 

where Fmn is the Maxwell tensor and F^i is dual to it. Consider now these master equations in the case of a static 
spherically symmetric spacetime. 

B. Static spherically symmetric non-minimally coupled fields. Reduced three-parameter system of 

equations: one-parameter models 

1. Preliminaries 

Using Schwarzschild coordinates, the line element for a static spherically symmetric system can be put in the form 

ds^ = B(r) c^dt^ - A{r) dr^ - {d0^ + sin^ 9 dip^) , (10) 

where B and A are metric potentials that depend on the radial coordinate r only. This form of the line element is 
useful as when B and 1/A are simultaneously zero it signals the presence of an event horizon. Assume also that the 
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electromagnetic field inherits the static and spherical symmetries. Then the potential four- vector of the electric field 
Ai has the form 



(11) 



From ifTTj) the Maxwell tensor is equal to Fik = ^o(?')(<5[(5^ — S^S"^), where a prime denotes the derivative with respect 
to r. To characterize the electric field, it is useful to introduce a new scalar quantity E{r) as E'^{r) = —^FmnP"^^- 
Then the electric field squared is E'^{r) = ^^(^o)^ from which one obtains in turn Fro = —{AB)i E{r). Since the 
expressions 1/A and ^/ AB enter frequently in the master equations, it is sometimes convenient to use the functions 
N{r) and a{r) defined as N{r) = and cr(r) = A{r)B{r) . In summary, the functions 



E\r) 



and 



N{r) 



A{r\ 



, a{T) = ^A{r)B{r) 



are alternatives to the functions ^o(f'); -^{f)-, and B{r). 



(12) 



(13) 



2. Key equation for the Maxwell field and its solution 
The Maxwell equations ^ with (fTTI) give only one non-trivial equation, namely, 

[r2i?(r)(l + 2x°V(r))]' = 0, 
which can be integrated immediately to give 

E{r) jr^ 1 + (qi + 92 + qz) (n" + 3iV' ^ + 2N^ 



(14) 



+2r(2<zi + (^N' + iV^^ + 2qi{N - 1)| = Q , (15) 

where Q is a constant, to be associated with the central electrical charge of the solution. This equation gives the 
electric field of a central charge, corrected by the radial component of the dielectric permeability tensor, 1 + 2^^Q^{r). 
This component describes the vacuum screening effect on the central charge, due to the interaction of the vacuum with 
curvature, analogously to the screening of a charge by a non-homogeneous medium in a spherical cavity. Supposing 
the spacetime to be asymptotically fiat, i.e., ^'^(^(oo) = 0, one can see that ifTSj) yields asymptotically the Coulomb 
law E ^ Q/f'^, and the constant Q indeed coincides with the total electric charge of the object. 



3. Key equations for the gravitational field 

The equations for the gravitational field Q with jSJ-llHl) and the metric potentials lfTO|) redefined as in (fT3|) can be 
rewritten as a pair of equations for N{r) and (j(r), respectively. 



[r{l-N)]' 



{E^y N{qi+q2 + q3) 



8N 



0(91 + 92 + <73) NV' + + — (2<Zi + 92 



N 



^E' 



r 

" N' 



N 



+ (gi + 92 + 93) N" + -iN' - + 2N 2- 



a 



r 



-(291 + 92) ( 2— + 2 — - + — 



r a 



91- 



(iV-1) 



(16) 
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2a' 



= ^{Ey{q,+cn + q,) + {E^)' 



(91 



92 + 93 j 



+E^ 



4\ ^2 

- + 291 + 92 - 
r I r 



2a' 



I ^ ^ \- 293 
(91 +92 +93)— 2" 



(17) 



The first equation can be reduced to an equation for E{r) and N{r), by extracting the term ^ from the second one. 
The second equation contains the unknown functions E(r) and a{r) only. Thus, Eqs. lfT5 |) - (fT7|) form the key system 
of equations for the non-minimal Einstein-Maxwell model of a static spherically symmetric object. It is a system of 
three ordinary differential equations of second order for the three unknown functions, E{r), N{r) and a{r). The form 
of equations is not canonical. In principle, the electric field E{r) can be extracted explicitly from (flS)) as a function 
of N" , N', N, ^ and r. Inserting such E{r) into the Eqs. (Il6])-(Il7l), we obtain equations for N{r) and a{r) of 
fourth order in their derivatives. 



4- General features and notes 

Below we focus on models admitting solutions to Eqs. (fT5| - lfT7|) . such that they can be represented by a series 
expansion regular at r — 0, i.e.. 

Air 0) = AiO) + A'{0) r + ^A"iO) + ... (18) 

where A{r) symbolizes generically the functions E{r), N{r), and a{r). Our purpose is to find solutions satisfying three 
conditions: First, the electric field E{r) should be a continuous function regular at r = {E{0) ^ 00) and also should 
be of Coulombian form at r ^ 00. Second, the metric functions N{r) and a'^{r^N(r) should be regular at r = 0. Third, 
in terms of the functions A{r) and B(r) the asymptotic flatness requires that B"{oo) = A" {00) = i?'(oo)A'(oo) = 0, 
and A{oo) = const, B{oo) = const. So, essentially, one can put, a{oo) — 1 and iV(oo) — 1. 

Note that the regularity of the functions E{r), N{r), and a{r) at r = does not guarantee that the solution of 
the Einstein-Maxwell model is characterized by regular curvature invariants. For instance, when N(0) is finite but 
N{0) ^ 1, the model displays a conical singularity and the scalar invariants of the curvature tend to infinity as r ^ 0. 
In considering solutions such that the fields are finite we try to be as close as possible to Bardeen's idea of having 
black hole solutions without singularities, by finding a regular E{r) and putting the metric coefficients in the form 
N{r) = 1 — 2Mr^{r^ + rQ)~i and a{r) = 1, for some ro [l3|- As we will see it will turn out that this is not achieved, 
since although the electric and metric potentials are regular, the black hole solutions found here are singular at the 
center, where the curvature invariants blow up. Notwithstanding, these solutions are very interesting. Using the 
ansatz (fT8| and the Eqs. lfT5|) - (fT7|) one can couple the values E{0), N{0), a{0), and 91, 92, 93. The relations are 
different for the cases a{0) = and cr(0) 7^ 0, which we now analyze. 

(i) cr(0) = 0; When a{0) = 0, but cr'(O) ^ 0, one obtains from the system (Il5l)-((l7|) that (^^j (r ^ 0) i, and the 
decomposition ifTSj) is valid at r — > 0, when the following conditions are satisfied, 

39192 + gi + 29193 p2/n^^ , ^ 1 n 

E[0) — ■ =Q, S (0)(9i +92) - 1 = 0, 

391 +92-93 

N{0) [2(391 + 92 - 93)] = 291 + 92 , (19) 

for generic 91, 92, and 93. There are two specific cases. When 91 + 92+93 = 0, but both 291 + 92 ^ and 91+92 ^ 0, 
then Q = ■^E{0){q2 — qi) and iV(0) = j. When 91 + 92 + 93 — and 291 + 92 — 0, simultaneously, then KqiE^{0) = — 1, 
providing 91 is negative, and iV(0) is fixed by the relation A^(0) = 1 + 2qi'E{o) ^- ^® ^^^^ '^(^) 7^ ^' below, 

here the Ricci scalar R{r) takes an infinite value at r = 0. 

(a) a{0) 7^ 0; When all three functions, E{r), N{r), and a{r) are regular at r = 0, and a{r), appearing in the 
denominator of Eqs. (fT5l) - ifT7|) . does not vanish at r = 0, one obtains from the system lfT5|) - (fT7l) the following set of 
equations 



S(0) 29i[iV(0)-l] = g, ^2(0) 293 = 0, 
7V(0) [1 + nE^iO)iqi - 92 - 293)] = 1 + 91 '^^^'(0) , 



(20) 
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for generic qi, q2, and q^. Since the charge of the object, Q, is considered to be non- vanishing, one obtains immediately 
from the first equation of the set that E{0) ^ and N{0) ^ 1. Thus, we infer, 

Q ^ 2[N{0)-l] + ^E{0)Q ^ 

2E{0)[N{0)-1]' 2N{0)kE^0) ' ?3 " ' ^ > 

The relations (|20l) give that the curvature invariants are infinite in the center r = 0. For instance, when N{r) and 
cr(r) are regular in the center and cr(0) ^ 0, then the Ricci scalar 



R{r) = 2N— + N" + 37V'— + - [ N' + N— + —{N - 1) (22) 
a (7 r \ a ) 

tends to infinity at r -> 0, since A^(0) 7^ 1, as well as generally 7V'(0) ^ -Ar(0)^(0). 



The order of differential equations and the choice of the parameters: one-parameter models 



Now we want analyze the simplest cases of the system of Eqs. (fT5| - lfT7|) . One sees that there is an immediate 
simplification when 91 + (?2 + 93 = 0, since second order derivatives and products of first order derivatives disappear 
from the equations. In such a case the system ifTSjl - lfTT]) reduces to 



E{r) 



2r{2qi + 92) [N' + N— 



2qi{N-l)\=Q 



(23) 



— =r(2(7i+92)(S')'-g3i^' 



(24) 



^E' 



Ml^..(i.^)'(2,,+,2)A^ 

K 

qi + 2r(2gi + (72) ( N' + N-] + N{3qi + 92) 



(25) 



Moreover, the three-dimensional matrix, composed of the coefficients in front of the first derivatives {E'^)' , N' and 
a' , has rank two. This means that for gi + (72 + 93 = the system l(23|l - (|25|l can be reduced to one algebraic equation 
and two differential equations of the first order. The corresponding algebraic equation is 



K(2qi + q2){r^ - 2qi)E^ ~ 2«;(2gi + 92)^^' + P + 2q:i{N ~ l)]E - Q = , 



(26) 

which in turn links two functions, E{r) and N{r). Now from equation l(26|) . one can consider three subcases, that 
emerge from the case gi + 52 + 93 = 0. First we consider briefiy the trivial case in this context, qi = q2 = q^ = 
q, q — 0, i.e., the Reissner-Nordstrom limit, see subsection IIIIBI Then we consider two interesting non-trivial cases: 
first, qs 7^ 0, second, 53 = 0, When q^ ^ it is easy to express N{r) in terms of E{r). In the subsection IIII CI 
we consider a specific model in this class, characterized by the supplementary condition qi = —q,q2 — 2q, q^ = 
—q (i.e., (?i + 52 + 93 = 0, 2qi +92 = 0) . This is the Gauss-Bonnet type model, which has been considered as an 
important model in [1, [^, and for which we present an extended analysis with relevant new details. For the second 
case qi = ~q, q2 = q, qs = (i.e., 91+92+93 = 0, 93 = 0) , E{r) decouples from N{r) and we deal with a cubic 
equation for the determination of the electric field. We will consider such a model, the integrable model, in subsection 

niDi 



III. SOLUTIONS OF THE REDUCED THREE-PARAMETER MODEL: SOLUTIONS OF 

ONE-PARAMETER MODELS 

A. Basic quantities and variables 

One should first define three quantities, tm, rq, Eq, as follows 



rM^2GM, rQ = VG\Q\, Eq = 



Q 



(27) 
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Now, the models we are going to discuss here are essentially one parametric, with qi, 52, and being a multiple of 
some parameter q. It is useful to introduce first a quantity Vq, given through 

Tg = ^2^, and 2q = ±rl, (28) 

with Tq being a radius. From tm, tq, and r^, one can then construct two independent dimensionless quantities, 
namely 

a=^ = ±4, (29) 

and 

if ^ ^ . (30) 

The a quantity gives the deviation from the standard Reissner-Nordstrom case, and K fixes the ratio between the 
total mass to the total charge of the object. 

In addition, for what follows below, it is useful to write the equations of motion by defining two dimensionless 
variables, a normalized radius x and a normalized electric field Z{x), defined as follows 

r E(r) 
x=-, Z{x)^-y. (31) 

Also, the function N{x) can be defined in terms of another useful function y{x), i.e., 

N{x) = 1 - ^ . (32) 

X 

The physical interpretation of y{x) is connected with the effective mass of the object, i.e., M(r), as we will see below. 
With these quantities defined we now discuss the Reissner-Nordstrom limit and the two new models. 

B. The Reissner-Nordstrom limit: qi = q2 ~ qi = q, with q — 

When the non-minimal parameters qi, q2, and q^ are set to zero, i.e., 

qi=q2=qz^q, q = 0, (33) 

and so a = as well, we can integrate immediately Eqs. (fT5|) - (fT7l) . In terms of the above functions, the system of key 
equations can be rewritten as 

x^Z-l = 0, (34) 
xN'{x)+Nix) = l-2Z + x^Z'^, (35) 
a'ix) 



The solutions to these equations are 



. (36) 



Zix)^^, (37) 



, , K ~l/x 
N{x) = 1 '— , (38) 



a{x) = 1. 



(39) 
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The function y{x) in (|32l) is here given by 



yix) ^K-^, (40) 



where K is defined in ((3Q|) . and one can see through Eq. l(35|) that y obeys the equation 

dy 1 



Of course one can transform to the original fields E{r), A{r), and B{r), giving 



E{r) = ^, (42) 



1 2GM{r) , , 

1 ^ , (43) 



A{r) 



2GM(r) , , 

B{r) = 1 ^ , (44) 



which are the usual Reissner-Nordstrom functions. Note, then, that the physical interpretation of y{x) given in l(40|) 
is connected with the effective mass of the object, i.e., M(r), which in turn is given by the definition 

Nir)'^l-'-^ = l-yM. (45) 
^ ' A{r) r X ^ ^ 

Thus, y{x) is related to the dimensionless effective mass, since y{x) = ^ , with M{r) = AI — M being the 
asymptotic mass of the object. 

We study in this context the usual static spacetimes with a — 0, i.e., the Schwarzschild and Reissner-Nordstrom 
spacetimes, which are special solutions of the full system of equations. These two solutions are heavily singular, both 
the metric and the Kretschmann scalar diverge at r = 0. So these solutions are outside the spirit of the solutions 
we want to find. They do not serve as models. However, they are of interest to set the nomenclature, and to have a 
gauge with which we can compare the solutions we find in the two models studied below. Note that for these solutions 
a^{x) = A{r)B{r) = 1. In these cases the problem of searching for horizons is equivalent to finding their radial 
position Th through the solutions of B{r) = = 1 — ^'^^^('') = q, where M{r) is the effective mass. In terms of 
dimensionless mass y{x), see Eq. (|45| . or Eq. l(40|) . this condition can be written as y{x) — x. For the Schwarzschild 
metric y{x) = K, where K = rM/rq is defined in ((30l) . Then, one obtains only one horizon at x — Xh = K, which is 
just the Schwarzschild radius, r ~ tm- For the Reissner-Nordstrom metric one has y{x) — K — ^, and the equation 
K — ^ = X gives three different cases: (i) K > 2 (i.e., tm > Srg, or GM^ > in the standard notation): there 

are two solutions x^ = ^ + \j ^ ^ 1; ^ind Xh = -f- — — 1, corresponding to the outer and inner horizons, 

respectively, of a usual Reissner-Nordstrom black hole, (ii) K — 2 (i.e., tm = 2rQ, or GAP = Q^): there is one 
horizon only, given by Xh = 1, corresponding to an extremal black hole, (iii) K < 2 (i.e., tm < Srg, or GM^ < Q^): 
there is no solution to the equation y{x) = x. The object is a naked singularity. 

C. The Gauss-Bonnet model: qi = — g, §2 = 2g, qa = —q (i.e., qi + q-z + qs = and 2qi + 92 = 0) 

1. Preliminaries 

Consider now, in Eqs. lU-lISl), the following specific one-parameter model 

qi = -q, q2 = 2q, = -q (i.e., qi + q2 + qs = 0, 2qi + q2 = 0) , (46) 

for some parameter q. In this model the susceptibility tensor is proportional to the double-dual Riemann tensor and 
is divergence-free i.e., Xikmn = Q *^ikmny ynX**^™" = 0- Moreover, the coupled Einstein and electromagnetic 
equations are second order in the derivatives, which is the reason why this model is called a Gauss-Bonnet model. 
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Gauss-Bonnet gravity in five and higher dimensions has the property that its equations are of second order, Lovelock 
gravity being a generaUzation of it. Actually, one can show that the model specified by (|46l) comes from Kaluza-Klein 
reduction to four dimensions from five a dimensional Gauss-Bonnet theory, i.e., Einstein gravity plus a Gauss-Bonnet 
term Q. Using l|46p . Eqs. l[23|) - (|25| convert, respectively, into 



E{T){r^ + 2q{l-N)] = Q, 



(47) 



r — = nqE 
a 



(48) 



(49) 



After appropriate redefinitions these equations agree with the ones discussed in The correspondingly modified 
algebraic Eq. l(26|) coincides with (|47|) . We now consider two ways of analyzing this system of equations: first, we use 
a power series expansion, and second, we apply the formalism of dynamical systems. 



2. Abel equation and its solution 



(I) The Abel equation 
Using from Eq. (gT]) that 



N{r) 



Q 

E{r) 



one obtains the Abel equation (see, e.g., [23]) for E{r) from (|49|) . 

rE'{r) = E- 3— -B^ - KqE^ . 



(50) 



(51) 



Similarly, eliminating E{r) from Eq. ([47)1 one can transform Eq. (|49|) into the Abel equation for a function 9(r), here 
defined as 



Thus, Eq. (|49| is given by 



e(r) = 1 - N{r) . 



[re'(r) + e(r)] [r^ + 2qe{r)] = 



(52) 



(53) 



Clearly, the values E{0) and N{0) are related through N{0) = 1 — 2qE(o) ■ Searching for solutions N{r) regular at 
r = 0, we have to consider E(0) to be non- vanishing, E{0) ^ 0. In such a case Eq. ((5T|) yields nqE^ifi) — 1. This 
means that q has to be positive in these solutions. Thus, although we analyze all cases, we tend to focus in models 
with a = 4f > 0. Finally, a possible solution, regular at r = 0, should be characterized by a-(O) — 0. It is convenient 

to use the auxiliary quantities Vq = v^Sjgj, rg, and a used before (see Eqs. l(28|l - (|29|) ). to write 



E{tS) = 



«9 TQi-q 2 \ q rq y/a 



(54) 



Using Eqs. l(47j) - (|49l) . plus the asymptotic conditions cr(cxD) = 1, N{oo) — 1, as well as the condition that the electric 
field is asymptotically Coulombian E Q/r^, one can obtain the following formula for the asymptotic mass M, 



TM = 2GM = lim 



r(l - N)~ -tkQE + 2NKqE^' 



= lim r(l - TV) , 



(55) 
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(II) The solution of the Abel equation for small r 

In the vicinity of the point r = the solutions for E{r)^ N{r) and o-(r) are assumed to have a polynomial form of the 
type given in Eq. ifTS]). The decomposition of a regular solution E{r) with non- vanishing E(0) is 



Q 



E{r) 

and the corresponding N(r) with finite N{0) is given by 

1 



(56) 



N{r) 



1 



-(-V 



For this solution the effective mass M{r) becomes equal to zero at r = 0. Moreover, -/V(O) 
non-minimal radii coincide, i.e., rg = r^, a= 1. 

(Ill) Power series expansion with respect to ^ 
The decomposition of the electric field yields 



(57) 

0, when electric and 



1-a 



(?)-!:*(?) 



n-l 



(58) 



where the 6„ are defined below. Infinity is a regular point for iV(r), thus, taking into account (|55l) one obtains the 
following decomposition of N{r) 

„2 



1 



A{r) 



n— 5 



4rr 



n+l 



where again, the 5„ are given below. The function cr(r) = A{r)B{r) is equal to one in the Schwarzschild and the 
Reissner-Nordstrom cases, but not in general. When g 7^ the logarithm of this function can be represented by the 
decomposition 



In 



^ ' 4 V r / ^ n+3 \ r J 



n=4 



00 00 ^ 



n— 4 m— 4 

The 6n coefficients can be taken from Eq. l(59|) . and starting from 65 can be found by the recurrence formula 

n-l 



with 



m=0 

-1 , 62 = 63 = 
&6 = , 67 



n + 3 



(60) 



(61) 



, 64 = a- — , 
4rQ 



a 
5 ' 



2_VM_ 

7r2 



bg = a' 



9rM 
32rQ 



(62) 



60)) are regular at r = cx) and absolutely converge in the interval r > tq H{a), where H{a) 



^ i. Note that the terms 60 and 61 are the Schwarzschild and Reissner-Nordstrom terms, respectively, and 



The decompositions ([58 
lim„^oo 

that the 5„ for n > 2 are the terms that give the post-Reissner-Nordstrom behavior. Note also that for g = 0, i.e., 
the Reissner-Nordstrom case (or the Schwarzschild case when, further, (5 = 0), the function a in l(60|) is equal to one, 
as it should. For r ^ 00, Eqs. l(58l) - l(60|) give useful asymptotic formulas, showing that, for E{r) and N{r), the first 
post-Reissner-Nordstrom terms are of fifth order in (-^), and that the decomposition for log cr starts with a term 
of fourth order. When necessary one should convert from N and a io A and B. Numerical calculations, see Figs. 
[DlSl confirm that the corresponding curves tend to the corresponding horizontal asymptotes, when r goes to infinity, 
i.e., E{r) tends to zero, and and B{r) tend to 1. It follows from ((SS)) that, for positive q, q > Q, the curvature 
coupling effects on the electric field are analogous to a dielectric medium, since the asymptotic electric field effectively 
decreases. For negative g, g < 0, there are no solutions with N{r) regular at r = 0, and, since we are mostly interested 
in regular or quasi-regular solutions, we do not fully discuss this case. 
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FIG. 1: Non-minimal solution of the Gauss-Bonnet model, with qi — —q, q2 — 2q and gs = — g, of gravitational electrically 
charged objects characterized by a — 2g/rQ, for K — 2^/2 [K > 2) - Plots (a), (b), (c), and (d) depict the metric potentials 
1/A{r) and B{r), the function log ^y A{r)B{r), and the electric field E{r)/EQ, respectively, as functions of a; = for solutions 
with difi^erent values of the non-minimal quantity a. The Reissner-Nordstrom black hole has a = 0, the curves of which are 
clearly shown in the plots. In this case the curves for 1/A{r) and B{r) have two zeros representing the inner and outer horizons, 
and for r ^ oo they go to one, while the curve E{r)/EQ tends to zero, respectively, and when r these curves tend to — oo 
(the electric field is in a logarithmic scale). For a < and < a < ao the black holes behave quite similarly as the case a = 0, 
with two horizons, and the function E{r)/EQ tends to finite values as r ^ 0. For a — ao the curve for 1/A{r) tends to a finite 
negative value when r — > 0, and takes the value zero only once. On the other hand B{r) has two zeros one at the same point as 
1/A{r), the event horizon, and the other at r = 0, signaling the presence of a singularity there. The function E{r)/EQ tends to 
finite values as r 0. For a > ao the curves and B{r) have one zero, and thus one horizon only, at the same r, and then 
tend to infinity as r ^ 0, an analogous behavior to the Schwarzschild black hole. The function E/Eq tends to finite values as 
r — > 0. See text for more details. 



3. The dynamical system associated with the model 
We now study this model, specified through Eq. (|46l) . using a dynamical system analysis. 
(I) First analysis: The plots and numerics 

(A) Key dynamic equation: In order to find the regular solutions E{r), N{r) and a{r) in the whole interval < r < oo 
let us transform the master Eqs. (|47l) . l(48|) and l(49|) to the independent variable dimensionless radius, and 

to the unknown dimensionless function y{x) given in Eqs. ((38|) and l(40|) . i.e.. 



y{x) = x[l ~ N{x)] 



(63) 
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FIG. 2: Non-minimal solution of the Gauss-Bonnet model, with qi — — g, q2 — 2q and gs = — g, of gravitational electrically 
charged objects characterized hy a — 2g/r-Q, for .K" = 2 - Plots (a), (b), (c), and (d) depict the metric potentials 1/A{r) and 
B{t), the function log ^ A(r)B(r), and the electric field E{r)/EQ, respectively, as functions of a; = for solutions with 
different values of the non-minimal quantity a. The extremal Reissner-Nordstrom black hole has a = 0, the curves of which are 
clearly shown in the plots. In this case the curves for 1/A{r) and B{r) have one double zero representing an extremal horizon, 
and for r ^ oo they go to one, while the curve E{r)/EQ tends to zero, respectively, and when r ^ these curves tend to — oo 
(the electric field is in a logarithmic scale). For a < and < a < ao the black holes behave quite similarly as the case a = 0, 
with one horizon, and the function E{r)/EQ tends to finite values as r — > 0. For a = ao the curve for 1/A{r) tends to a finite 
negative value when r — + 0, and takes the value zero only once. B{r) has two zeros one at the same point as 1/A{r), signaling 
there is only one horizon, and the other at r = 0, signaling the presence of a singularity there. The function E{r)/EQ tends to 
finite values as r ^ 0. For a > ao the curves -^p^ and B{r) have at the same r, one zero, and thus one horizon only, and then 
tend to infinity as r ^ 0, an analogous behavior to the Schwarzschild black hole. The function E(r)/EQ tends to finite values 
as r ^ 0. See text for more details. 



The physical interpretation of y[x) is connected with the so-called effective mass of the object, M{r), see Eq. (|45| . 
Putting these definitions into Eq. (|49|) . we obtain the following key equation 

dyix) X 

= TT • 64 

ax x'^ + ay[x) 

This equation is indeed a key one, since using its solution, y{x), we can represent explicitly the electric field by 

Eix)^E^Zix), with Zix).-^^±^^, (65) 
the metric function N{x) by Eq. I|45p . and <j{x) by the integral form 

\iia{x)^\ii^A{x)B{x) = a / —Z^{x'). (66) 

J CO 
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FIG. 3: Non-minimal solution of the Gauss-Bonnet model, with qi — —q, q2 — 2q and = — g, of gravitational electrically 
charged objects characterized hy a — 2q/r^, for K = 1 {K < 2) - Plots (a), (b), (c), and (d) depict the metric potentials 1/A{r) 
and B{r), the function log ^J A{r)B{r), and the electric field E{r)/EQ, respectively, as functions oi x = for solutions with 
different values of the non-minimal quantity a. The Reissner-Nordstrom naked singularity has a — Q, the curves of which are 
clearly shown in the plots. In this case the curves for 1/A{r) and B{r) have no zeros, and for r — > oo they go to one, while 
the curve E(r)/EQ tends to zero, respectively, and when r — > these curves tend to —ao (the electric field is in a logarithmic 
scale). For a < and < a < ao the naked singularity behaves quite similarly as the case a = 0, and the function and E{r)/EQ 
tends to finite values as r ^ 0. For a = ao the curve for 1/A{r) tends to a finite positive value as r — > 0, while B{r) has a zero 
at r = 0, signaling the presence of a singularity there. The function E{r)/ Eq tends to finite values as r 0. For a > ao the 
curves and B{r) have one zero at the same r, and thus one horizon only, and then tend to infinity as r ^ 0, an analogous 
behavior to the Schwarzschild black hole. Thus, by tuning the non-minimal quantity a one can turn a Reissner-Nordstrom 
naked singularity, which has a = 0, into a black hole, when a > ao. The function E{r)/ Eq tends to finite values as r — > 0. See 
text for more details. 



Moreover, A follows from A — 1/N and B from B = cr'^ /A, i.e., 



B{x) 



1 



exp|2a y""^z2(a;')| 



X 

Note that the function y{x) is also a function of the quantity o, so in general should be written as o). 



(67) 



(B) Three typical cases: In this problem there are two independent dimensionless quantities, constructed from rj\/, 
rg, and r^, namely a and if, see Eqs. l(29|) and ((30l) . Note that K, besides fixing the ratio between the total mass 
and the charge of the object, gives the value of the dimensionless mass y{x) at r = oo, since y{oo) = K . Taking into 
account the quantity K, and in conformity with the Reissner-Nordstrom solution, let us distinguish three different 
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FIG. 4: Reduced mass profiles t/(0, a) and y{x, a) for the non-minimal Gauss-Bonnet model, with q\ = — g, §2 = 2g and ga = — g, 
of gravitational electrically charged objects characterized by a = 2g/rQ - Plots of the curves i/(0, o) and y{x,a) are shown, 
with a; = r/rg. For ao < a (with ao > 0) the characteristic curve y = —x^/a is not intercepted by the function y{x,a). For 
a = ao the intersection takes place at a; = 0. For < a < ao the point of crossing floats along the characteristic curve, and the 
integral curve has two branches. Finally, for a < there is no intersection of y{x,a) with the characteristic curve. In the plot 
2/(0, a) it is shown explicitly the existence of a point obeying y(0, ao) = with positive ao. 



situations, (I) K > 2, (ii) K — 2, (iii) K <2, within each situation the quantity a can vary from zero to infinity. Figs. 
[His] display typical cases in each situation, and Fig. |4] shows the behavior of y{x). In slightly more detail: (i) K > 2 
(i.e., VM > 2rQ): For K > 2,-we use K — 2\/2 as a typical value for the numerical analysis, see the plots in Fig. [U (see 
Fig. [IK for Fig. [TJd for B{r), Fig. for In ^y A{r)B{r), and Fig. [iji for E{r)). When a = this case gives the 

usual Reissner-Nordstrom black hole with two horizons. For other as there are also black holes, some with different 
properties, (ii) K = 2 (i.e., tm = 2rQ): For if = 2 see the plots in Fig. [21 (see Fig. [2k for Fig. [2b for B{r), 

Fig. for In ^/A{r)B{rj, and Fig. [H for E{r)). When a = this case gives the extreme Reissner-Nordstrom black 
hole with one horizon. For other as there are also interesting solutions with black holes, (iii) K <2 (i.e., tm < 2rQ): 
For K < 2, we use X = 1 as a typical value for the numerical analysis, see the plots in Fig. [31 (see Fig. [3K for -jjyj, 

Fig. [3Jd for B{r), Fig. [3b for In ^yl(r)B(r), and Fig. [Sji for E{r)). When a = this case gives a Reissner-Nordstrom 
naked singularity, a solution without horizons. For other as there are also solutions. 

(C) Scaling of the key equation: The key Eq. fM)) remains invariant after the following scale transformations: 

x^^x, y^Ky, a->-^a. (68) 

Thus, the critical values of a that one may eventually encounter when K — 1, are also critical values that one can 
easily find for arbitrary K using the formula ao{K) = ao(l). 
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(II) Second analysis: Critical properties of the family of the solutions 

(A) About the mass function y{x,a) when x = 0, y{0,a) - the critical value of the quantity a, oq; We now write 
explicitly that y is a function of both x and a, y — y{x, a) since this is important to our analysis. Plots of this 
dimensionless mass function y{x,a) are displayed in Fig. [H In Fig. UK, y{x,a) is shown for several values of a, 
and in Fig. |4b, a plot for y{0,a) as a function of a is shown. A simple qualitative analysis shows that the mass 
function y{x, a) at the central point x = 0, i.e. ?/(0, a), as a function of the quantity a, has to possess a zero. Indeed, 
when a = 0, y(x,Q) = K — ^, (where, recall, K = 7^), corresponding thus to the Reissner-Nordstrom solution. 
When a — 00, y{x, 00) = K, corresponding thus to the Schwarzschild solution. Note as well that when a ~ —00, 
?/(x, —00) = constant. In addition, when x = 00, ?/(oo, a) = K, & condition at infinity that holds for arbitrary a. In 
order to prove our assertion, that ?/(0, a) as a function of a possesses a zero, consider then y(0, a) as a function of the 
quantity a in the interval < a < 00. One can see, that y(0, 0) — —00, and y{0, 00) — K > 0. Supposing that y{0, a) 
is continuous in such an interval, one can conclude that there exists at least one specific value of the a quantity for 
which 2/(0, ao) — 0, where ao is the value of a for which y(0, ao) — 0. Fig. IHd shows that, for X = 1 the zero of the 
function y(0, a) = happens, when a = ag ~ 7.49. But the most interesting fact is that the curve y = t/(0, a) displays 
a discontinuity in the first derivative with respect to a just at a = oq. One can see explicitly a finite jump of the 
derivative at this point, a — qq. Nevertheless, the function y — y{Q, a) itself is continuous at this point. That is why 
Oq is a critical value of the quantity a. For other Ks the critical values can be found using the scaling properties (|68l) . 
yielding ao{K) — For instance ao(2) ~ 0.468 and, for the typical case we study, ao(2v^) = 0.117. For negative 
a, 2/(0, a) ^ asymptotically when a —00. 

(B) The critical points of the associated autonomous two-dimensional dynamical system: The key Eq. (|64l) can be 
put as an autonomous dynamical system 

y — x, x — x^ + ay, (69) 

where ' = ^ and r is an auxiliary parameter. In Eq. l(69|) there is one critical point at 

(a;,y) = (0,0). (70) 

In the vicinity of this critical point the variables x and y are connected by the relation x'^ — ay"^ = constant, which 
means this point is a saddle point when a is positive, and a center when a is negative. If a > there are two 
separatrices y = ±:^. The equation for N{r) — -^^^ can also be written as dynamical. 



iV = (1 - TV) + 2(7(1 -N f 

which is much more complicated than Eq. (|69l) for y{x). 
the critical points immediately, 

(r, N) = 



-r|, f = r[r2 + 2q(l-7V)] , (71) 
Nevertheless, if q is positive (i.e., a is positive), one can find 

^0,1±^). (72) 



In order to present the integral curves for the total interval of the auxiliary parameter r, we resort to numerical 
calculations. The results are presented in Figs. [lEl It is clear that for the critical a = aq the curve for N{r) — -^y^ 

tends to one when r 00, and takes a finite value ^^^y = 1 — at the center of the object, r = 0. This critical 
curve is a separatrix between the curves having a > ao and those having a < gq. The same type of behavior happens 
with the curves for B{r). One also has that for a = ao there exists a unique integral curve for y{x) with asymptotic 
value given by y{oo) = and for r — > one has y{0,ao) = 0. Again, this curve behaves as a separatrix. Similar 

reasoning goes to the curve E{r). In other words, the critical point N = 1 — given in Eq. (|72|) . is just a saddle 
point, corresponding to the critical value ao, obtained numerically, and it is the final point of the unique integral 
curve, which coincides with the separatrix of the function y — at small x. When q < and, thus, a is negative, 

there are no regular or quasi-regular solutions to the master equations for the whole interval < r < 00. 

(C) Vertical asymptotes and the electric barrier: From Eq. (|64|) one sees that, when y{x,a) — — the derivative 
y'{x, a) becomes infinite and vertical asymptotes appear in the graph of y{x) versus x. At the point x = x* , for which 
y{x* ,a) — — the electric field l(65l) becomes infinite, and so these vertical asymptotes can be interpreted in terms 
of an electric barrier. 
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When a < 0, Figs. [TJi,[2li, and[3li show that for all three values of K the curve for the electric field has a form of 
a finite barrier. Indeed for r ^ cx) the electric field E{r) tends to zero, then at some radius it reaches a maximum, 
the barrier height, and finally goes to a finite positive value when r ^ 0. A charged test particle with sufficiently 
high energy can overcome this barrier, be trapped in the potential well and then oscillate inside. When a — 0, more 
specifically, when a ^ 0_, the barrier height increases and tends towards the center r ^ 0. For a — (see again Figs. 
[Hi, [2li, and [S]!) one obtains the standard, Reissner-Nordstrom, behavior, E{r) — Sj. This means that the electric 
barrier has become infinite, taking its maximum value (an infinite value) at the center r = 0. In other words, the 
vertical asymptote for the electric field appears at r = 0. Finally, when a > this vertical asymptote and the position 
of the infinite electric barrier shift from the center towards positive r values, then stop at some value of the quantity 
a, drift back again to smaller values of r and, finally, the infinite electric barrier disappears at a > oq, (with oq > 0). 
When a — ao the curve tends to the horizontal asymptote at r ^ 0, there is yet no trap. When a > ao, one can see a 
finite electric barrier with the corresponding traps, for all values of K. 

Since the integrals in l(66|) and l(67|) diverge when E(r) is discontinuous, then in searching for solutions with regular 
functions E{r), -^pj, and -B(r), we should reject all the cases where vertical asymptotes appear. For a > ao, numerical 
calculations show that vertical asymptotes do not appear. So we will mainly consider solutions for this range of the 
quantity a, i.e., a > ag. 

(D) Horizons: To analyze the a > ag case, with a non-singular electric field and possible non-singular metric potentials 
1/A and B, we have studied previously, for comparison, the usual static spacetimes with a = 0, i.e., the Schwarzschild 
and Reissner-Nordstrom spacetimes, which are special solutions of Eq. (fTO|) . One has for these that A{r) B{r) = 1. 
For a > ao one sees from l(60|) that A{r) B{r) ^ 1, in contrast to the a = case. Nevertheless, as in the a = case, 
horizons are still given by the condition -^pj = ov y{x) — x. We have analyzed this numerically. For a > ao the 

results are the following: (\) K > 2 (i.e., tm > Srg)): A typical case is K = 2\/2, see Fig. [H One finds that for 
a > ao = 0.117, there is only one horizon, (ii) K = 2 (i.e., rM — 2rQ): For K — 2, see Fig. [21 One finds that for 
a > ao = 0.468, there is one horizon also, (iii) K < 2 (i.e., tm < 2rQ)): A typical case is K — 1, see Fig. [Sj One 
finds that for a > ao = 7.49, there is only one horizon also. 

Note that Figs. [l][3|show that when o > ao, for all formal possibilities {K > 2, K = 2, K < 2) the curves 1/A tend 
monotonically to minus infinity and cross the line 1/A = only once. Thus, for arbitrary a > ao the plots of 1/A are 
continuous, irregular at the center and characterized by one horizon. These solutions have thus an analogous behavior 
to the Schwarzschild black hole. Moreover, by tuning the non-minimal quantity a one can turn a Reissner-Nordstrom 
naked singularity, with K < 2 and a = 0, into a black hole, when a > ao. Now, since the a = ao is a very special case, 
we discuss it in particular. 

(E) The solution with a — ao: In the framework of the model under discussion, i.e., when qi — — g, q2 = 2q, 
and <73 = —q, all the three functions, E{r), and B{r), are regular in the interval < < oo if and only if 

a — ao^K) = (recall a = ^). This means that we deal with a one-parameter family of exact solutions, the 

arbitrary quantity being K — and all the curvature coupling constants, qi, q2 and (73, being expressed explicitly 

via K. The critical q corresponding to ao is thus, with the help of Eqs. (l27|)-(l30l), given by, qag = 7.49Q'^/32GM"*. We 
now consider these solutions in more detail: (i) There are two different solutions, corresponding to the separatrices 
y{x) — at small x. The physical solution, the solution that gives the appropriate limit when a; 00 and has no 

jumps on the derivative of the characteristic curve ay{x) = —x^, is y{x) ~ +"^^- The plots of ^4^, B{r) and E{r) 
are displayed in Figs. [HSj (ii) The solutions, which we discuss, are characterized by finite values at the center, and 
E{0) ^ 0, A{0) ^ 0, and A^(0) ^ 0, but a{0) = and B{0) = 0. (iii) For this model ^ B{r), and there are two 

distinct critical radii; first, the radius for which — N{r) — 0; second, the radius for which B{r) = which signals 
an infinite redshift surface, and an event horizon in the case of static spacetimes, such as the ones we are treating 
here, (iv) For — 0, one finds that such a solution exists when < rg, i.e., ao < 1; this is a necessary condition 

(see, e.g., fSTj) ). Moreover, within this case it is possible to have such a zero when K'^ > 7.49, i.e., K > 1.65 or 
TM > 1.65rQ. For Vq > rg, one has, = 1 — ^ is positive and without zeros, (v) For the infinite redshift surface 

and event horizon, B{r) — 0, one finds that when a = ao, for the critical quantity, the function B{r) is zero both at 
r = and at the radius for which = 0. (vi) The Kretschmann scalar diverges at r = 0, so although the metric 
functions are regular, spacetime is not, time stops. 

(F) Remark: Some important aspects derived from qualitative and numerical analyses of this non-minimal model, with 
qi = —q, q2 = 2q and 93 = —q, can be found in Q. Our results are in thorough concordance with this initial analysis. 
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We have supplemented those aspects on several grounds, of which we stress briefly three novel details obtained here: 
(i) We have found a complete converging decomposition of the solution of the Abel equation based on the recurrence 
formula This gives us not only the asymptotic decompositions for r ^ oo, but also the possibility to link the 

limiting formulas for r ^ and r ^ oo (see l(58|) - ([60|) ). (ii) We have formulated and discussed the problem of the 
inflnite electric barrier, associated with the vertical^asymptote appearing when + ay{x, a) = 0, thus completing 
physically and mathematically the analysis given in [;!?]. (iii) We have found, first, the critical value of the non-minimal 
parameter g, namely, = 7.49(5^/32G'Af'*, second, the scaling law of the critical parameter for different values of 
the asymptotic quantity K = tm/tq, namely, qa„{K) = 7.49/if , third, the significance of the choice for K, namely, 
K > 2, K — 2, and K <2,m the qualitative analysis. 

D. The integrable model: q\ = —q, q2 ~ qz ~ Q (i.e., Qi + 92 + 53 = and ga = 0) 

1. Preliminaries 

Now, we consider the model with 

qi = -q, q2 = q, qs = (i.e., qi + q2 + qs = 0, q^ = 0) , (73) 

which, as we have seen, has the property that E{r) decouples from N{r) and we deal with a cubic equation for 
the determination of the electric field. Since the basic feature of the model is that it is integrable, we call it the 
integrable model. For such a model the susceptibility tensor can be written in terms of the Einstein tensor 

Qik ^ pnk_ijigik foiio^j^g^ ^^fcm« ^ | \p^^ gkn _ Qin gkm _^ Qkn gtm _ G<fe™g»n] ^ ^^ere we have put qi = -q2 = -q. 

Thus, the model becomes one-parametric and we can introduce the dimensionfull quantities defined above, rg, Eq, 
Tg, and the dimensionless quantities a and K. We assume that Eq inherits the sign of the charge Q and the quantity 
a can be positive or negative depending on the sign of q. Then, we introduce the two dimensionless variables, the 
normalized radius x and the normalized electric field Z{x), defined in l|3ip . In terms of these, the system of key 
equations can be rewritten as 

a{x^ + a)Z^ -2aZ^ - x'^Z +1 = 0, (74) 
xN'{x) +N[1- ax(Z'^y{x)] =l~2Z+{x'^+ a)Z^ , (75) 

^^-a(z7(x). (76) 

Clearly, Eq. l(74|) for Z{x) is the key equation for finding N{x) and a{x) from l(75|) and (|76l) . respectively. If instead of 
N{x) one uses y{x) then Eq. l(75|) gives an equation for ^ of the type of Eq. I|4ip or Eq. l(64|) . but more complicated, 
which for this analysis is not very illuminating. For this model it is better to start analyzing the electric field E{r), 
or its redefinition Z{x). 

2. Electric field 

Consider now Eq. l(74|) in detail. Eq. ((74l) is a one-parameter algebraic equation of third order for the dimensionless 
electric field. Below we denote its solution as Z{x,a). The function Z{x,a) can be generally presented by the 
well-known Cardano formula, nevertheless we prefer to analyze qualitatively this one-parameter family of solutions. 
Depending on the value of the quantity a the solution Z{x, a) can possess one or three real branches. The corresponding 
plots are presented in Fig. [H When a — 0, one obtains, as it should, the Coulombian solution Z{x,0) — The 
curves displaying Z{x, a) for non- vanishing values of the quantity a are more sophisticated. 

(i) a > 0; When a is positive, the functions Z{x, a) take finite values for all values of the quantity a, see the curves 
a, b, c, d in Fig. [H The initial values Z{0,a) satisfy the cubic equation a'^Z^{0,a) — 2aZ'^{0,a) -1-1 = 0. There is 
only one real solution of this cubic equation for a > 0, if the discriminant V = j^j^{27a — 32) is positive, i.e., when 
II < a < oo, see box a of Fig. \E\fov details. This plot displays three real branches of the solution Z{x, a), nevertheless, 
only one of them is defined on the whole interval < x < oo. Two other branches are defined for x > Xnunio^) only 
and contact at the point x = x^in{a). Only one branch tends to the horizontal asymptote Z = at a; ^ oo. When 
< a < II , the discriminant V is negative or equal to zero, which guarantees that there are three real starting points. 
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FIG. 5: It is shown, in the non-minimal integrable model, with qi = —52 = ~q and 93 = 0, the rescaled electric field Z{x,a), 
{Z{x, a) = E{r)/ Eq), as a function of the scaled radius x {x = r/rg) of gravitational electrically charged objects characterized 
by a = - (a) Displays the solution Z{x,a) when the non- minimal quantity a satisfies the inequality || < a < 00. There 

are three real branches of the solution Z{x,a). Nevertheless, only one of them is defined on the whole interval < x < 00. 
Only one branch tends to the horizontal asymptote Z — at x ~* co. (b) Displays the solution Z{x, a) in the case 1 < a < ||. 
There are three real starting points Zi{0,a), ^2(0, a), ^3(0, a) for the three corresponding branches of the solution Z{x,a). 
Nevertheless, two branches are not continuous, only the third being defined on the whole interval < a; < ca. (c) Displays 
the solution Z{x,a) for the important case a = ao = 1, the Fibonacci soliton. (d) Displays the solution Z{x,a) in the case 
< a < 1. (e) Displays the solution Z{x,a) in the case a = which is a Coulombian electric field. The curve is not continuous, 
(f) Gives an example of the solution for negative a with a vertical asymptote. When a tends to zero remaining negative, the 
vertical asymptote shifts towards the line a; = 0, and the solution Z{x, a —> 0-) converts, finally, into the Coulombian solution. 
At large values of x the plots of the function Z{x, a < 0) tend to the Coulombian curve Z{x, 0) for all values of the quantity a. 
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Zi{0, a), ^2(0, a), ^3(0, a) for the three corresponding branches of the solution Z{x, a), see the curves in boxes b, c, d 
of Fig. m Nevertheless, when 1 < a < two branches of the solution Z{x, a) are not continuous, only the third being 
defined on the whole interval < a; < 00, see box b of Fig. [5l When < a < 1 all three branches are continuous and 
are defined on the whole interval < 2: < 00, one of them is asymptotically Coulombian. There are three horizontal 
lines Z = i=, Z — and Z — +— i==, which yield distinct ranges for the functions Zi(x, a), Z2(x, a), and Z^ix, a). 

Clearly, the curve of Coulombian type is in-between the separatrices Z = and Z = +— The model with critical 

VH 

a, call it ao again, is the one that has a = 1, i.e., a = ao — l.This model can be solved analytically, and we consider 
this case in detail below. Finally, when a tends to zero remaining positive, the starting points Zi{0,a) tend to minus 
infinity, and ^2(0, a), .^3(0, a) grow infinitely. Clearly, at a ^ 0+ the branch ^2(2^, a) is the only branch that remains 
visible at finite values of x, and is Coulombian. 

(a) a — 0; When a = 0, one obtains the Coulombian solution Z{x, 0) — Jj-, see the curve in box e of Fig. [H 

(Hi) a < 0; When a is negative, the coefficient (a;^ + a) in the first term of Eq. I|74p vanishes at a: = y^joj and the 
line X = ^/\a\ is the vertical asymptote of the graph Z{x,a), see the example of the curve for a = —1 in box f of Fig. 
m When a tends to zero remaining negative, the vertical asymptote shifts towards the line x = 0, and the solution 
Z{x, a ^ 0_) converts, finally, into the Coulombian solution. At large values of x the plots of the function Z{x, a < 0) 
tend to the Coulombian curve Z{x, 0) for all values of the quantity a. Thus, for a < the solutions Z{x, a) are not 
regular at all in the range < a; < 00. 



3. Metric functions 



To study the gravitational part of the solution, we observe that the solution of l|75l) for N can be presented in 
quadratures 



N{x,a) = 1 - ie''^'(^''') {K- 



ie+a)Z^ta) - 2Z(e,a)+2a^Z(e,a)-Z(e,a) 



(77) 



where Z{x, a) is supposed to be already found. The constant of integration, K, can clearly be related to the asymptotic 
mass of the object M. When a; ^ 00 and Z{x 00, a) ^ Eq. ([77)) yields that K = as defined in l(30|) . On 
the other hand, searching for a solution N{x, a), which is finite at a; = 0, and taking into account that Z{0, a) is finite, 
we should require that 



K 



dC 



In this case the formula l|77|) transforms into 

iV(a;,a) = -e'^^'("''^) T d^e-'^^'(«''^) [{^^+a)Z^{^,a)-2Z{C,a)+l] , 



providing 

7V(0, a) = 1 + aZ^{0, a) - 2Z(0, a) . 
Since the electric field at the center satisfies the condition 

1 - 2aZ^{0, a) + a^Z^iO, a) = , 



(78) 



(79) 



(80) 



(81) 



(see l(74l) 1. then N{0,a) can be rewritten as N{0,a) = 1 — ^^^^ in accordance with the first relation from (|20|) . 
Thus, for the family of solutions with regular functions N{r, a) and E{r, a) the quantity if is a function of the value 
Z(0, a), i.e., depends on the quantity a according to the formula (|78|) . K — K{a). 
Also, from (|76l) for a one finds 



cr(x, a) — (t(oo, a) exp {—a [Z'^{x, a) — Z^(oo, a)] } , 



(82) 
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FIG. 6: Non-minimal solution of the integrable model, with qi = —q, q2 = q and qa = 0, of g ravi tational electrically charged 
objects characterized by a = -S?- - Plots (a), (b), (c), and (d) depict the functions A, B, log V AB and E/Eq, respectively, as 

functions of a; = for three typical values of the quantity a, when < a < 1. These functions are regular and take finite 

values at r = 0. The solution with a = ao = 1 is a solution without horizons, since and B are positive everywhere. Indeed 
the a = ao = 1 solution is a soliton of the model, the Fibonacci soliton. It has a mild conical singularity at the center, and is a 
solution without horizons. 



with Z{x,a) being found from ([74l) . For the asymptotically Coulombian branch we have to set Z(cx),a)=0 and thus 
tT(cx), a) = 1. Then, the function 1/A is taken from 1/A = N and the function B is taken from A and a, B = jA. 

In the Reissner-Nordstrom solution and in the previous discussed model, i.e., the Gauss-Bonnet model (see Sections 
IIII Bl and IIII C| respectively), we divided the solutions according to K > 2, K — 2, and K <2. Here it is no more 
convenient to make such a division. The reason is that K is not a free quantity here, rather K — K{a). Since it is 
not a free quantity, we cannot classify the models with respect to it, we can only calculate this quantity (numerically) 
after solving the problem as a whole. In the Gauss-Bonnet model of section llll CI we could classify through K because 
the equations for Z{x) and N(x) cannot be decoupled, thus, Z(0, a) and A^(0, a) are connected. This means that 
we can choose if as a convenient quantity for the classification, the dependent quantity being Z(0, o). Here, in this 
model, the solution for Z{x) satisfies a decoupled equation, the latter does not depend on N{x). Thus solving the 
decoupled equation for Z{x, a), we can classify the quantity Z{0, a) as the independent one. Then, we obtain N{x, a) 
and, as we see, -/V(0, a) depends on E(0, a), an so if is a dependent quantity, K — K[a). An explicit example where 
K{a) is calculated is given below for the case a = — 1. So, as when discussing the electric field, we again divide 
the analysis into three cases, here with subcases. 

(i) a > 0; We should further divide into three subcases: 

oo > a > 1: The electric field is discontinuous or irregular at the center. Since in this work we focus on regular electric 
fields everywhere, although interesting, we do not discuss these models here. 
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a = 1, i.e, a = ao = 1 : We plot in Fig. [6l boxes (a), (b), (c), (d), the functions 1/A, B, log V AB, and E/Eq, 
respectively. From the figure it is clear that a=l is characterized by the absence of horizons. For a = oq = 1 the 
solution is regular, or better, quasi-regular, and is a soliton of the theory, the Fibonacci soliton. Due to its interest, 
the case a = ao = 1 will be solved next explicitly. 

1 > a > 0: This case is interesting. In this case we plot in Fig. \E[ boxes (a), (b), (c), (d), the functions l/A, B, 
log V AB, and E/Eq, respectively, for two positive values of a within this range, namely, o = 0.999,0.750. From the 
figure it is clear that positive a in this range gives one horizon, and the solution possesses a quasi-regular center. 
Since the singularity at the center is a conical one, these black hole solutions can be considered as quasi-regular 
solutions, and thus are of great interest. Note that extremal black holes have only one zero, which in turn is a double 
zero. So from the figure above, the non-minimal black holes do not characterize as extremal. Rather, they are of the 
Schwarzschild type, with one horizon, a spacelike singularity, but here, different from Schwarzschild, the singularity 
is mild, it is a conical singularity. 

(a) a = 0; It is the Reissner-Nordstrom case. For this case the electric field is irregular at the center. 

(in) a < 0; For these cases the electric field is discontinuous or irregular at the center, like the Reissner-Nordstrom 
case, and we do not discuss these models here. 



4- Exact solution: a = ao = 1, the Fibonacci soliton 



Let a = ao = 1. Then the cubic equation l(74|) takes the form 

{Z - 1) [(1 + x'^)Z^ + {x^ - 1)Z - 1] = . (83) 

One sees that Eq. (|83)) splits into one linear equation and one quadratic equation. One branch of solutions 
of l(83l). the linear one, describes a constant electric field Zconst(a;,l) = 1, or, equivalently, E{r) = Eq. This 
branch is of no great interest. Another branch ^nonCouiomb(2^, 1) is given by the function ^nonCouiomb(2:, 1) = 
2(^i^rci^ [l — — Vx"^ + 2x2 ^ 5j ^ -which is bounded. The graph of this function starts from ^nonCouiomb(0, 1) = 

— ^2~^ tends asymptotically to the line Z = —1. The behavior of such electric field is not of Coulombian type, 
and will be not discussed further. Yet, there is a third branch. The branch ^couiomb(a;, 1) = Z{x,l) given by the 
formula 

^(^' 1) = onl fl - + + 2x^ + sl , (84) 
2(1 -I- x^) L J 

describes a Coulombian type electric field. At x ^ oo, one has Z{x, 1) jj-, or equivalently, E{r) — > ^. The graph 
of this function starts from Z{0, 1) — ■ Interesting to note that the starting points Z(0, 1) are associated to the 
well-known Fibonacci series and the "golden section" = ^^'^ = = 1.618.... For the Coulombian type solution, 

the function N{x, 1), which we write simply as N{x) when suitable, is regular at the center only if the constant K 
satisfies ([78l) . The corresponding quadrature for N{x) is 

Nix) = £ aiO [e + 3- , (85) 

where cr(x) is given by dM]). Clearly, N{oo,l) = iV(oo) = 1 and A^(0, 1) = iV(0) = so that l-iV(O) = i = 0-1, 

and the relations (|20l) are satisfied. The plot of the function N{x) for a = ao = 1 is shown in Fig. [6l Clearly, the 
function N{x) is positive in the interval < x < oo. For the Coulombian type solution l(84|) the function a{x,l), 
which we write simply as a{x) when suitable, is given by 



J 3 + il-x^Wxi + 2x^ + 5 + x^ 1 
= exp I 1 , (86) 

with (t(0, 1) = cr(0) being equal to exp{-(l + (j))}. Then one finds 1/A from 1/A = N and B = a'^/A. The function 
B{x) is also positive, and B{0) = (1 - i)e-2(i+'« ~ 0.002032. Thus, this solution is a solution without horizons and is 
regular. Although the curvature scalars diverge, the singularity at the center is a mild one, it is a conical singularity. 
The asymptotic mass M of the object defined as 

tm = 2GM = lim {r [1 - B{r)]} , (87) 
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is represented in this case by the integral 



Numerical calculations give the value 



which yields in addition, 



1 ^^'iO 1 



M ~ 0.442 



^(0 {e+3-v^^T2e+5) 



\Q\ 



K 



2VGM 



G 



0.884. 



(89) 



(90) 



This is thus a very interesting solution. It is a soliton, in the sense that is made of the very own fields of the theory, 
the gravitational and electric fields, it is a solution without horizons, and it is quasi-regular, with a conical singularity 
at the center. 

Note also that the value a = uq — 1 can be regarded as a critical one. There are two reasons for this. First, it is 
clear from Fig. \E\ that for the unique case a = qq = I there exists a bifurcation point, in which the two branches of the 
curve Z{x,a), or E{r,a), intersect. When a > 1, the Coulombian branch of the electric field curve is discontinuous. 
When < a < 1, three continuous regular branches exist. The second reason is that aX a = qq = 1 the corresponding 
curve on Fig. [6] plays a role of a separatrix; when < a < 1, desirable curves, continuous and regular at the center, 
exist, otherwise they do not appear. 



E. Summary of the results 

In Fig. [7] we present a table in which the results of the various models studied are summarized. 





qi=-q. ^2=2g, q^=-q 


qi = -q, ^2 = ^3 = 


K = \, 
= 7.49 


K = 2, 
= 0.468 


K = 242, 
=0.117 


i: = 0.884, 


a <0 


Finite electric field Eif ) 
Infinite l/A(0), B(0) 


Finite E(o), 
Infinite electric barrier 


Horizon | 1 Horizon | 2 Horizons 


a = 


Minimal Reissner - Nordstrom model 


<a 


Discontinuous electric field Eif ) 


Finite E{r) , llA{r) , B{r) 

1 Horizon 


a = 


Fimte E{r), l/A{r), B{r) 


Horizons 


1 Horizon 


1 Horizon 


Horizons 


a> 


Finite E{r) , infinite l/^(0) and B{0) 


Discontinuous E{t) 


1 Horizon 


1 Horizon 


1 Horizon 



FIG. 7: In this table it is diplayed the main results on the studied models. 



IV. CONCLUSIONS 



We have shown that the original non- minimal Einstein-Maxwell theory with three parameters gi, 52, and 53, is 
reducible in natural different ways to a theory with one parameter q only, in which the three parameters obey two 
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relations between themselves. We have then studied two special models for static spherically symmetric solutions 
obeying the following requirements: the electric field E{r) is regular everywhere in the interval < r < oo, being 
Coulombian far from the center. From the solutions of this class we extract the ones, for which the metric coefficients 
-^yj and B{r) are regular at the center r = and tend to one asymptotically as r ^ oo. 

The first non-minimal model, the Gauss-Bonnet model (with qi = —q, 92 = Sg, (73 — —q, q free), displays charged 
black hole solutions with one horizon only, when the dimensionless non-minimal quantity a, with a — "iq/rq naturally 
appearing in the model, exceeds a critical value ao, a > uq. Although the black hole is electrically charged the 
solutions found have one horizon only, and are similar in this connection to the Schwarzschild solution. When a < uq, 
the solutions are discontinuous in the interval < r < cx3, or irregular at the center r — 0. Another main result in this 
model is that there exists a unique solution, the solution for the critical value a=ao, which does not possess horizons 
and is characterized by regular fields E{r), 1/A{r) and B{r), with B{0)=0 and A{0) ^ 1, although the curvature 
invariants blow at the origin. 

The second model, the integrable model (with qi = —q, q2 = q, qs = 0, q free), is also characterized by one critical 
value a = ao = 1 of the non-minimal quantity a. When a < or a > 1, the solutions are irregular. When < a < 1 one 
obtains black holes with electric field regular everywhere and with only one horizon, like the Schwarzschild solution. 
Finally, when a = oq = 1, i.e., at the critical value of the quantity a, there exists a solitonic solution with a conical 
singularity at the center, but otherwise well behaved. This solution can be called the Fibonacci soliton, since the 
well-known cj) number {cj}— — 1.618), associated with the golden section, appears naturally in the expressions for 
the central values of the electric field G\Q\E{0)=(j> and the metric coefficients are also related to cj), namely 1~^^=^ 

andB(O)=(l-i)exp[-2(l+0)]. 

Summing up, we can say that the non- minimal curvature induced interaction between the gravitational and electro- 
magnetic fields provides an electric field, of static spherically symmetric charged objects, which is regular everywhere 
for different relations between the coupling constants qi , 52 and 93 . As for additional regularity of the metric co- 
efficients B{r) and A{r), the non-minimal interaction can provide models which have very specific, critical, values 
for the coupling constants, in which the geometry has at most a conical, and thus mild, singularity. This is in line 
with problem posed by Bardeen [l7|, where one should look for theories with regular black hole solutions. We have 
partially solved it within these models. 
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